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The conversion of multilayer graphenes into sp 3 -bonded carbon films on metal surfaces (through 
hydrogenation or fluorination of the outer surface of the top graphene layer) is indicated through 
first-principles computations. The main driving force for this conversion is the hybridization between sp 3 
orbitals and metal surface d z 2 orbitals. The induced electronic gap states and spin moments in the carbon 
layers are confined in a region within 0.5 nm of the metal surface. Whether the conversion occurs depend on 
the fraction of hydrogenated (fluorinated) C atoms at the outer surface and on the number of stacked 
graphene layers. In the analysis of the Eliashberg spectral functions for the sp 3 carbon films on a metal 
surface that is diamagnetic, the strong covalent metal-sp 3 carbon bonds induce soft phonon modes that 
predominantly contribute to large electron-phonon couplings, suggesting the possibility of 
phonon-mediated superconductivity. Our computational results suggest a route to experimental realization 
of large-area ultrathin sp 3 -bonded carbon films on metal surfaces. 



Diamond is known for its extraordinary thermal conductivity and mechanical hardness, and the production 
of large-area crystals has been a long term goal. The conventional synthetic methods employed for the 
production of artificial diamond usually require extremely high pressure and temperature 1,2 , or plasma- 
enhanced chemical vapor deposition (CVD) 3 " 5 . Advances in the growth of large-area graphene 6,7 and bilayer 
graphene through multilayer graphenes 8 " 10 suggest another possibility: the chemical conversion of a few layers of 
graphene into sp 3 -bonded carbon films. For example, previous theoretical results showed that the complete 
hydrogenation 11 or fluorination 12 of a single-layer of graphene can yield a thermodynamically stable sp 3 -bonded 
carbon layer. Indeed, partial hydrogenation of monolayer graphene 13,14 and evidence of complete fluorination of 
monolayer graphene membranes 15 have been experimentally achieved. Similarly, the conversion of a few layers of 
graphene has attracted attention, and hydrogenation of the outer surfaces of bilayer graphene 16,17 to achieve all 
sp 3 -bonded carbon films has been treated theoretically. We suggest however that such an approach will be 
challenging experimentally as extensive hydrogenation or fluorination was stated to be needed on both the top 
and bottom surfaces. 

Here we present a promising alternative, namely the conversion of graphene layers on a metal surface into sp 3 
carbon films through surface hydrogenation or fluorination. From the practical point of view, the commonly used 
methods for growing graphene films on a metal substrate, either through CVD or as a result of the transfer of 
grown graphenes 6 " 10,18 , are already providing appropriate target systems for the study of the sp 2 to sp 3 conversion 
of carbon films on metal surfaces. Our computation results suggest a method of preparing an ultrathin 'diamond' 
layer over even very large area by conversion of appropriate multilayer graphene on metal surfaces. The resulting 
structure could also be viewed as a carbon-based metal-insulator junction that may reveal novel two-dimensional 
(2D) phenomena, as we suggest in this study, or as a new type of electronic material 19 . 

Results 

We first investigated the effect of hydrogenation or fluorination on the free surface of the AB-stacked (Bernal-type 
stacking) bilayer graphene. Simulating experimental generation of gas phase atoms 13,20 , hydrogen or fluorine 
atoms were considered to be available on the surface. Figure la is a schematic of the bilayer graphene before the 
adsorption of such gas atoms (denoted by I). Experimental studies reported that the hydrogen chemisorption 
coverage of a single-layer graphene is less than half of the carbon atoms 20 . It was however suggested that bilayer or 
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Figure 1 | The energetics related to the adsorption-induced geometric 
changes of bilayer graphene. The geometry (a) before (group I) and (b) 
after the hydrogen adsorption (group II), and (c) after formation of 
interlayer C-C bonds (group III), (d) Relative energies for each 
configuration and the energy barrier in unit of kcal mol -1 of hydrogen 
adatoms. The insets in (d) show three different atomic arrangements of 
hydrogen adsorption on the outer surface of graphene (boat-like, zigzag- 
like, and chair-like) in the group II. Only the chair-like configuration is 
available in the group III with an energy of about 23.5 kcal mol -1 higher 
than that of the boat-like configuration in II. TS represents the transition 
state. The larger gray and smaller blue balls represent the carbon and 
hydrogen atoms, respectively. 

multilayer graphene can have more extensive hydrogenation, which 
can be up to half the coverage of the outer surfaces 16,21 . In Fig. 1, our 
super-cell model contains hydrogen (fluorine) atoms that could cover 
half the carbon atoms on the outer surface. After the chemisorption 
of hydrogen (fluorine) atoms, the two graphene layers can be sepa- 
rated by the Van der Waals' (vdW) distance, as shown in Fig. lb, or 
can form interlayer C-C bonds, as shown in Fig. lc. The overall 
energetics for the cases of hydrogen adsorption is described in Fig. Id. 

The half coverage hydrogenation of the surface of graphene can 
exist in various configurations. Top views of the patterns, named as 
'boat-like', 'zigzag-like', and 'chair-like' 1112 , are presented in the 
insets of Fig. Id. In the boat-like configuration, two hydrogen atoms 
are located at the two ends of the bridge of the C-C bond. In the 
zigzag-like configuration, each carbon atom in every second zigzag 
chain is bonded to a hydrogen atom. In the chair-like configuration, 
the C-H bonds are equally spaced. These three configurations, in 
which the carbon layers are separated by the vdW distance, are 
together categorized as configuration II in Fig. 1. Only the chair-like 
configuration induces all-carbon sp 3 hybridization in the carbon 
atoms, leading to sp 3 dangling bonds on the inner side of the top 
graphene layer. This drives formation of covalent C-C bonds 
between the upper and lower graphene layers, as shown in Fig. lc 
(denoted as III). In contrast, dangling bonds are essentially not pre- 
sent for the boat-like and zigzag-like configurations owing to the 
pairing of two neighboring p orbitals in the form of local n bonding 22 . 
Thus, the boat-like and zigzag-like configurations are much more 
stable than the chair-like one either in the configuration II or III. 
Overall, the energetics shown in Fig. 1 indicate that the interlayer 
C-C bonds are not favored. Previously, it was shown that 



hydrogenation of only one side of graphene can cause the graphene 
to roll-up 23 . Such a structural instability of a flat graphene upon one- 
side hydrogenation further reduces the chance for the formation of 
interlayer C-C bonds. 

When the one-side hydrogenated graphene bilayer is positioned on 
transition metal surfaces, the energetics for conversion to sp 3 films are 
significantly different than for the isolated bilayer case. In Fig. 2, we 
consider the same adsorption-induced conversion of bilayer gra- 
phene on the (0001) surface of Co (a few other metals are also 
considered, as discussed below). Besides the existence of Co below 
the graphene layers, the adsorption configurations of C-H bonds 
presented in Figs. 2a, b, and c, and the insets of Fig. 2d (configura- 
tions I, II, and III) are the same as those shown in Fig. 1. The favored 
adsorption site of carbon atoms is found to be atop Co atoms at the 
interface, as shown in Fig. 2a, which agrees with previous studies 24 ' 25 . 
The formation energy differences between the boat-like, zigzag-like, 
and chair-like configurations in group II are nearly the same as those 
presented in Fig. 1. However, it is noteworthy that the presence of Co 
greatly stabilizes the configuration III; this is in obvious contrast to the 
situation without the metal substrate. This stabilization occurs 
because of the saturation of the otherwise unstable sp 3 dangling 
bonds with Co surface states. Further, as shown in Fig. 2d, the trans- 
ition barrier from configuration II to configuration III is negligibly 
small (about 0.5 kcal mol" 1 ). 

Bond lengths and bond angles in configuration III reveal obvious 
sp 3 carbon (diamond-like) bond features. The calculated interlayer 
C-C and intralayer C-C bond lengths are 1.58 A and 1.52 A, respect- 
ively, and these values are comparable to the sp 3 bond length (1.54 A) 
in bulk diamond. The tetrahedral angles of the C-C-C bond within 
the layer and that formed with the interlayer C-C bond are 110.2° 
and 108.5°, respectively; these values are close to the tetrahedral bond 
angle of 109.5°. Fluorination of the graphene bilayer on Co(0001) 
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Figure 2 | The energetics of conversion of bilayer graphene into sp 3 
carbon film on metal surface, (a-d) The same energetics and geometric 
changes as the ones that were shown in Fig. 1 but with Co(0001) 
underneath the graphene bilayer. In contrast to the cases of Fig. 1, the 
configuration III (the chair-like adsorption) is now more stable than the 
boat-like configuration in II by an energy difference of about 8.1 kcal 
mol -1 . The pink balls are the cobalt atoms. 
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revealed similar features as those observed for the hydrogenated 
cases. The details of structural properties are provided in Supple- 
mentary Table SI in the Supporting Information. 

To understand the underlying electronic origin of the sp 2 to sp 3 
conversion, we investigated the electronic structures of the converted 
sp 3 carbon films on Co (0001). The projected electronic band struc- 
ture and the projected density of states (PDOS) of carbon atoms in 
the outermost layer are presented in Fig. 3a. Pristine single-layer 
graphene has the cone-like energy bands touching the Fermi level 
at the K point of the 2D Brillouin zone. Such characteristics, origin- 
ating from the delocalized n electronic states, largely disappeared 
upon the conversion from sp 2 to sp 3 bonding. The top of the valence 
bands, in which the p x and p y orbital states are degenerate, and the 
bottom of the conduction band, which is characteristic of a non- 
degenerate p z orbital state, are set far apart from each other, resulting 
in a direct gap of about 3.18 eV at the T point. The first-principles 
calculations yielded similar wide band gaps for graphane 1112 and 
diamane 16 ' 26 (which are sp 3 carbon structures with full hydrogenation 
of top and bottom surfaces of graphene and AB- stacked bilayer 
graphene, respectively). The minority spin states, presented in 
Fig. 3b, reveal the same features; thus, there is no measurable mag- 
netization in the outermost layer. 

Figure 3c shows the band structure and PDOS in the majority spin 
state of the carbon layer at the interface with Co. The corresponding 
results for the minority spin states are shown in Fig. 3d. The in-plane 
p bands depicted in orange (p x ) and green (p y ) are almost identical to 



the corresponding projected bands of the outer layer (Figs. 3a and b). 
However, the hybridization features between the carbon p z and 
cobalt d z 2 orbitals are prominent throughout the energy level, leading 
to the drastic changes in energy and dispersion of p z bands. For 
comparison, the projected electronic structures of majority and 
minority spins in the interface Co atoms are presented in Figs. 3e 
and f, respectively. Such strong hybridization causes exchange-split- 
ting in the carbon layer. The magnetic moment in the first carbon 
layer is calculated to be about —0.08 |i B per carbon atom and anti- 
parallel to the Co spin direction, which is twice the value observed in 
experiments for graphene adsorbed onto the Co(0001) surface 27 . 
These metallic states developed at the interface can penetrate into 
the sp 3 carbon layers near the interface, constituting the metal- 
induced gap state (MIGS), as discussed in the next paragraph. 

To estimate the range of the MIGS, we considered the possibility of 
converting graphene multilayers of various thickness (see further 
discussion on thickness-dependent energetics below) to sp 3 carbon 
films; for example, we considered the case of five graphene layers 
stacked in the AB configuration on Co(0001), as shown in Supple- 
mentary Fig. SI a. The PDOS resolved into each layer (see Supple- 
mentary Fig. Sic) shows that the first carbon layer (LI) at the inter- 
face with Co has electronic states persisting over the Fermi level with 
large exchange-splitting. The evanescent tails of such gap states decay 
rapidly and almost disappear in the third layer (L3). This indicates 
that such sp 3 C films with a thickness of only approximately 0.5 nm 
from the metal surface can likely provide good electric isolation. 
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Figure 3 | Electronic structures of the two-layer sp 3 carbon film on Co(0001). (a) Majority and (b) minority spin band structures and PDOS of carbon 
atoms on the top layer. The same (c) majority and (d) minority states projected onto carbon atoms at the interface. The same (e) majority and (f) minority 
states projected onto the top Co layer at the interface. The Fermi level is set to zero energy. 
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We next investigated the effect of the coverage of hydrogen on the 
transformation of AB-stacked graphene bilayer into an sp 3 carbon 
film. Figure 4a shows the formation energy (solid circles) and inter- 
layer distance (open circles) between the two graphene layers on 
Co(0001) as a function of the hydrogen coverage x of the outer 
(i.e., top) surface. The magnitude of the formation energy increases 
gradually with the coverage. A similar trend was observed previously 
in diamane structures 16 . The interlayer distance changes dramatically 
at around x = 1/3, as the interlayer C— C bonds develop at a higher 
concentration, whereas the two carbon layers remain at the vdW 
distance at a lower concentration. In Fig. 4a, to provide a generic 
model for the effect of the distribution of hydrogen adatoms, we 
assumed that the hydrogen adsorbates are distributed uniformly over 
the graphene surface. In realistic situations, there can be various 
other configurations with a given concentration. An example is pro- 
vided in Fig. S2 in the Supporting Information. With the same con- 
centration of x = 1/6, hydrogen adsorbates can be aggregated in a 
narrow region (Fig. S2a) or scattered uniformly (Fig. S2b). In the 
aggregated region (x = 1), the hydrogen adsorbates induce the C-C 
interlayer (sp 3 ) bonds in that region. This is consistent with the right 
side (x > 1/3, highly concentrated region) of Fig. 4a. Various adsorp- 
tion configurations, within the same adatom concentration, can lead 
to various total energies, as demonstrated in Supplementary Fig. S3. 
However, the effect of local distribution of hydrogen adsorbates can 
be described by the model given in Fig. 4a. When the coverage of the 
outer surface is less than one quarter of the maximal coverage (x = 
0.25 or C 8 H), the two carbon layers are separated by 3.12 A, which is 
essentially the interlayer vdW distance. For H coverage greater than 
half of the maximal coverage (x = 0.5 or C 4 H), the formation of the 
interlayer C-C covalent bonds is favored, and the carbon atoms on 
both layers predominantly have the sp 3 characteristics. This depend- 
ence on the coverage is largely consistent with Angus and Hayman's 
formulation 28 for the composition of sp 3 and sp 2 bonds in hydrogen- 
treated diamond-like hydrocarbons, as discussed in detail in the 
Supplementary Information. 

We now report the effect of thickness on the stability of fully 
hydrogenated (C 2 H) sp 3 carbon films on metal substrates. We con- 
sider two different stacking configurations of graphene layers, 
namely ABC and AA. The ABC-stacked graphene was converted 



into a cubic diamond-like layer, whereas the AA- stacked graphene 
transformed into a hexagonal diamond-like layer. The side views of 
the cubic and hexagonal diamond-like overlayers are presented in 
the insets of Fig. 4b: the cubic diamond consists of three alternating 
layers (denoted by a, p, and y) and the hexagonal diamond consists 
of two alternating layers (denoted by a and P). Figure 4b shows the 
formation energies of cubic and hexagonal diamond-like ultrathin 
films on the Co(0001) substrate. The negative formation energies 
indicate that the sp 3 -bonded diamond-like structures on the Co sur- 
face are thermodynamically and structurally stable for thicknesses of 
up to seven carbon layers. 

The preferable conversion energetics from sp 2 layers into sp 3 car- 
bon films can be applied to the cases of other metal surfaces. The 
reactivity of metal surfaces can be roughly classified into two groups, 
depending on the strength of their interaction with graphene: the 
physisorption group' (e.g., Cu, Pt, Al, Ag, or Au), corresponding to 
the surfaces on which graphene weakly physisorbs through vdW 
attraction, and the 'chemisorption group' (e.g., Ni, Co, or Pd), cor- 
responding to the surfaces on which graphene develops relatively 
strong chemical bonds through a hybridization between its p orbitals 
and the metal d orbitals 24 . In the present study, we selected Co(0001), 
Ni(l 1 1), and Cu(l 1 1) because of the convenience associated with the 
periodic super-cell calculation; the lattice mismatches of these metals 
with the graphene are only 1.6, 1.2, and 3.6%, respectively. As 
expected, the conversion energetics for graphene layers on Ni(lll) 
revealed features that were almost similar to those of Co(0001). 

A more remarkable point is the stabilization of the configuration 
III on the Cu(lll) surface. As discussed in Ref. 24, the overlay er of 
pure graphene (configuration I) and the one- side hydrogenated gra- 
phenes (configurations II) physisorb onto Cu(l 1 1) through the vdW 
interaction. However, for configuration III, the sp 3 dangling bonds at 
the interface develop covalent bonds with surface Cu atoms, similar 
to that observed for Co and Ni. The calculated relative energies of the 
hydrogenated bilayer graphene on Ni(lll) and Cu(lll) substrates 
(along with the case of Co(0001) and the case without metal sub- 
strate) are summarized in Table 1. Clearly, the stabilization by 
hybridization between the sp 3 dangling orbital and the d z metallic 
orbital in configuration III applies not only to the chemisorbing metal 
surfaces, such as Co and Ni, but also to the relatively less reactive 
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x in C 2 H X Number of carbon layers 

Figure 4 | The dependence of stability of sp 3 carbon films on the adsorption coverage and thickness, (a) The hydrogen coverage-dependent formation 
energy (solid circles) and optimized C-C interlayer distance (open circles) of bilayer graphene on Co(0001 ). The vertical dashed line represents the critical 
coverage (x ~ 1/3) beyond which the interlayer C-C bonds are predominantly favored. The insets show the optimized structures corresponding to the 
coverage just below (x = 0.25) and above (x = 0.5) the critical coverage. In this model, hydrogen adsorbates are distributed uniformly over the two- 
dimensional supercell. The definition of the formation energy is described in the Supplementary Information, (b) The thickness-dependent formation 
energies of the fully hydrogenated (C 2 H at the outer surface) sp 3 carbon films on Co(0001). Solid and open symbols refer to the formation of cubic and 
hexagonal diamond configuration, respectively. The insets show the side views of the optimized seven-layer diamond-like films on Co. The carbon atoms 
in different stacking layers (a(3y for cubic and oc(3 for hexagonal diamond configurations) are differentiated with gray scale. For simplicity, only the top- 
most Co layer is shown. Other atomic symbols follow the same convention used in Fig. 2. 
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Table 1 Relative energies of configurations 1, II, and 
each case with different metal, the energy reference 

Groups Configurations 


III depicted in Figs, 
s the configuration 

Pristine 


2a, b, and c on Co(0001 ), Ni(l 1 1 ), and Cu(l 1 1 ), respectively. In 
1. The energies are given in kcal mol" 1 of hydrogen adatoms. 

Cobalt Nickel Copper 


1 


0.0 


0.0 


0.0 


0.0 


boat-like 


-36.4 


-37.3 


-37.1 


-38.1 


II zigzag-like 


-35.4 


-36.5 


-36.2 


-37.2 


chair-like 


-15.6 


-16.2 


-16.4 


-18.8 


III chair-like 


-12.9 


-45.4 


-47.8 


-45.2 



metals, such as Cu. Very recently, Rajasekaran et al 29 reported 
experimental evidence for the partial conversion of a few layer gra- 
phene to sp 3 -bonded carbon layers on Pt(lll) substrate. 

The electronic structures of the sp 3 -bonded carbon layers on metal 
surfaces shown above revealed similar features as those of hole- 
doped diamond. This motivated us to study the features of elec- 
tron-phonon coupling that may lead to phonon-mediated supercon- 
ductivity. The superconductivity of hole- doped diamond (the 
substitution of B for C) has attracted broad interest 30,31 , and a possible 
modulation of critical temperature (T c ) in a single sp 3 carbon layer 
(i.e., graphane) through B doping has also been studied recently 32 . To 
investigate the electron-phonon coupling characteristics in the 
metal/diamond interface structures, we selected a two -layer sp 3 car- 
bon structure on Cu(lll) that possesses a perfect diamagnetic elec- 
tronic structure. The phonon DOS (PHDOS) of the hydrogenated 
sp 3 -bonded carbon layers on Cu(lll) is shown in Fig. 5a. For com- 
parison, the PHDOS of the two-layer sp 3 structure hydrogenated on 
both sides (i.e., diamane) is also presented in Fig. 5a. Both systems 
show quite similar PHDOS characteristics in the frequency modes 
ranging from 400 to 1300 cm" 1 ; the overall patterns are analogous to 
those of pristine diamond 33 . The sharp peaks around 1200- 
1300 cm" 1 can mainly be attributed to the shear motion between 
the carbon and hydrogen layer in the outer surface. The C-C in-plane 
and out-of-plane vibrations dominate the medium-frequency region, 
contributing to peaks at around 600, 800, and 1000 cm" 1 . Details on 
the optical phonon modes are given in Supplementary Fig. S4 in the 
Supporting Information. For the cases of B-doped graphane and dia- 
monds, the softening of these carbon stretching modes with increasing 
doping concentration has been discussed in relation to the electron- 
phonon coupling strength that is a key component determining T c 32,34 . 
By comparison with the diamane, a PHDOS characteristic of the sp 3 
layer on Cu is the presence of soft modes in the frequency range below 
250 cm" 1 that are due to the vibration of Cu-C bonds: in-plane modes 
with a peak at about 140 cm" 1 and out-of-plane modes with a peak at 
about 230 cm" 1 (indicated by arrows). 

Figure 5b shows the Eliashberg's spectral function a 2 F(oo) and 
the electron-phonon coupling constant calculated for the sp 3 - 
bonded carbon layers on Cu(lll). Here, X(co) is defined as 

X((o) = 2 f a ^-^- dco', and thus the total electron-phonon coupling 
J « 

o 

constant is X (co — > °°). The overall PHDOS characteristics dominate 
the a 2 F(oo) throughout the range of phonon frequencies. The calcu- 
lated X((6) increases in two steps at low and medium frequencies, at 
which the aforementioned Cu-C and C-C modes exert a critical 
effect. More specifically, a 2 F(oo) possesses a sizeable spectrum (indi- 
cated by an arrow) at the low frequency region, that extends up to 
250 cm" 1 , which can be attributed to soft phonons of the metal/sp 3 
carbon interface, which contribute considerably to the integrated X 
of 0.25. The high-frequency peaks in ot 2 F(co) at around 1200- 
1300 cm" 1 , originated from the shear motion of the outermost C 
and H layers, provides a negligible contribution to X 32 , because the 
integrand is the inverse of co. 

As a complementary study, in order to compare with the B-doped 
diamond systems, we also considered the effect of hole-doping on X. 



As is summarized in Supplementary Table S2, the quantity X 
increases with increasing B doping concentration and reaches 0.81 
at the doping concentration of 12.5%, which is substantially greater 
than the value of 0.56 obtained for a 10% B-doped diamond 34 . This 
sharp increase in X is mainly attributed to the enhancement of the 
electronic DOS at the Fermi level (N(E F )). These results imply that the 
X value can be engineered by choosing the metal substrate that can 
promote N(E F ) to the interface sp 3 carbon layers. To demonstrate this 
conjecture, we repeated our calculations for the sp 3 -bonded carbon 
layers on the (111) surface of Pt: the relatively large N(E F ) indeed 
leads to a large X of 0.45 without doping (see Supplementary Table 
S2). 

The computation for a thicker sp 3 carbon layer on metal surface is 
quite demanding, but we suggest that the results of electron-phonon 
calculations presented above should be valid for thicker diamond 
layers because the a -bonded Fermi level states in the sp 3 framework 
are rather localized in the vicinity of the metal/diamond interface. To 
verify the accuracy of the present computation, we calculated the 



a 0.08 




0 200 400 600 800 1000 1200 1400 
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Figure 5 | Phonon and electron-phonon coupling features of the sp 3 
carbon/Cu(lll) interface, (a) PHDOS and (b) Eliashberg spectral 
function a 2 F(co) and the electron-phonon coupling constant X((d) for the 
two-layer sp 3 carbon film on Cu(lll). In (a), the PHDOS for diamane is 
shown with the green line for comparison. The arrows indicate the C-Cu 
vibration modes. 
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electron-phonon features of known materials (i.e., MgB 2 ) and the 
results with references are shown in Table S2 in the Supporting 
Information. For the case of the sp 2 graphene overlayer on 
Cu(lll), X is quite small, and thus the possibility of phonon- 
mediated superconductivity is negligible. This is almost the same 
situation as the pristine single-layer and bilayer graphenes, as dis- 
cussed in the previous first-principles study 35 . Even though the cal- 
culated X and T c of our sp 3 carbon layers on metals are lower than 
those of MgB 2 , our results indicate that metal/sp 3 carbon interfaces 
are a new candidate for superconductivity studies. 

Discussion 

Through the first-principles density functional theory (DFT) calcu- 
lations, we investigated the energetics of transformation of a few 
layers of graphene into sp 3 -bonded thin overlayers on metallic sur- 
faces (Co, Ni, and Cu). Strong hybridization between the sp 3 dangling 
bonds orbitals and the metallic surface d z 2 orbitals stabilizes the sp 3 - 
bonded carbon layers. The dependence of the energetics of conver- 
sion on the surface coverage of C-H or C-F bonds and on the 
thickness of the stacked graphene layers was tabulated. Phonon spec- 
tra were calculated, and the features of electron-phonon coupling 
with the possibility of superconductivity were investigated through 
analyses of the Eliashberg's spectral function. Our results suggest that 
scaled preparations of multilayer graphene on metal substrates can 
lead to a fabrication route of very large area ultrathin sp 3 -bonded 
carbon films that would represent an entirely new carbon material. 

Methods 

The DFT calculations were performed using the projector augmented wave pseudo- 
potential method 36 , as implemented in the Vienna Ab initio Simulation Package 37 . 
Exchange and correlation interactions between electrons were described with the 
generalized gradient approximation formulated by Perdew, Burke, and Ernzerhof 38 . 
The long-range dispersion corrections for the interlayer interaction were taken into 
account within the semi- empirical DFT-D2 approach suggested by Grimm 39 . Spin 
polarization was taken into account for all the calculations. Some more details for 
supercell configurations, relaxation scheme, and /c-point samplings are given in the 
Supplementary Information. 

The PHDOS, Eliashberg spectral function a 2 F(co), and electron-phonon coupling 
^(co) were calculated using the linear response theory within the density functional 
perturbation theory 40 , as implemented in the Quantum ESPRESSO 41 . The critical 
temperature T c was estimated using McMillan's formula for the Eliashberg equation 42 , 



COin 

" 1.2 



exp 



-1.04(1 + A,) 
A,(l-0.62n*)-n* 



(1) 



where coi n is the logarithmic average frequency and u* is the screened Coulomb 
pseudo-potential. The values of T c shown in Supplementary Table S2 in Supporting 
Information were calculated with the typical choice of u* = 0.1 34 . 
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